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We present FLUKA and MARS simulation studies of the pion production and energy deposition 
in the Neutrino Factory baseline target station, which consists of a 4MW proton beam interacting 
with a liquid mercury jet target within a 20 T solenoidal magnetic field. We show that a substantial 
increase in the shielding is needed to protect the superconducting coils from too much energy 
deposition. Investigations reveal that it is possible to reduce the magnetic field in the solenoid 
capture system without adversely affecting the pion production efficiency. We show estimates of 
the amount of concrete shielding that will be required to protect the environment from the high 
radiation doses generated by the target station facility. We also present yield and energy deposition 
results for alternative targets: gallium liquid jet, tungsten powder jet and solid tungsten bars. 

PACS numbers: 13.20.Cz, 13.75.Cs, 14.60.Ef, 14.60.Pq, 29.20.-c, 29.25.-t 



I. INTRODUCTION 

The current baseline option for the Neutrino Fac- 
tory [1] is to use a 4MW proton beam interacting with 
a free-flowing mercury jet to create an intense muon 
beam The MERIT experiment has shown a proof- 
of-principle demonstration of a high intensity liquid mer- 
cury jet target [1]. The interaction of the proton beam 
with the mercury jet creates low-energy pions that are 
captured by the high field (~ 20 T) solenoid and trans- 
ported into a decay channel where the muons resulting 
from pion decay are collected and stored until they decay 
to neutrinos. 

In this paper, we present a series of simulation studies, 
using the FLUKA [j| and MARS @ computer packages, 
on particle production and energy deposition (radiation 
dose) calculations. We first show the simulation results 
for the so-called Study 2a geometry configuration, and 
then describe how the geometry needs to be modified 
to address various safety issues. We also compare use- 
ful muon yields and energy deposition doses for different 
target material alternatives. Finally, we present a study 
on the concrete shielding requirements that will be neces- 
sary to protect the environment from the high radiation 
doses emanating from the target station. 
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II. SIMULATION PARAMETERS 



As a starting point, the Neutrino Factory mercury jet 
target station geometry is based on the Study 2a con- 
figuration [6|] , as shown in Fig. [U with the appropriate 
20 T field map based on the dimensions and currents in 
the normal conducting copper (~ 6 T) and superconduct- 
ing coils (~ 14 T). Variations to the geometry are made 
to reduce the energy deposition in the superconducting 
magnets, as well as incorporating more engineering con- 
siderations to the overall target station design. In all 
simulations, the proton beam has a transverse Gaussian 
profile with a root mean square radius of 1.2 mm. The 
kinetic energy of the proton beam is nominally set to 
8 GeV, but is varied when finding the optimal number of 
useful muons from the target. The mercury jet is mod- 
elled as a simple cylinder with a radius of 4 mm, tilted 
at approximately lOOmrad to the magnetic z axis. We 
also investigate the yields and energy deposition for al- 
ternative targets, also tilted at lOOmrad to the z axis: 
liquid gallium jet, powder tungsten jet (50% density) 
and solid tungsten bars. The first two alternatives are 
also modelled as cylinders with a radius of 4 mm, while 
the solid target is modelled as a 20 cm long, 2 cm di- 
ameter cylinder with a density of 19.25 g/cc. The angle 
between the target and the proton beam at their inter- 
section (z — —37.5 cm) varies between 20 and SOmrad, 
depending on the initial kinetic energy, in order to opti- 
mise pion production Q. 



2 



150 
100 
50 



I 



-50 
-100 
-150 







SCI 


|--||— 1 SC2-SCS 


Air 

^ SC6-SC13 




Cu magnet 
WC ShieU 


T- 


- beampipi ' 
Fe plug 






steel casing 

Be window Vacuum 




\l{gpooj 















-500 



500 1000 

Z (cm) 



1500 



2000 



FIG. 1. Schematic of the Study 2a geometry of the Neu- 
trino Factory target system. The shielding is comprised of 
tungsten-carbide (80%) with water coohng (20%). The su- 
perconducting magnets are labelled SCn, where n = 1 to 13. 



III. STUDY 2A GEOMETRY 
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FIG. 2. The accepted muon yield per proton per GeV as a 
function of proton beam kinetic energy. The different points 
show what effect different simulation code versions have on 
the yield calculation results. ST2 and ST2a denote slightly 
different field tapering parameters. 



A. Muon yields 

One important figure of merit concerning the perfor- 
mance of the target is the total number of muons that 
pass through the Neutrino Factory cooling channel. This 
is calculated by first counting the number of pions, kaons 
and muons that are directly produced by the mercury jet- 
proton beam interaction. These secondary particles are 
then tracked through the solenoidal target decay channel 
up to z = 50 m, where they enter the muon cooling sys- 
tem. The ICOOL simulation package [8] is then used to 
find what fraction of these secondary particles end up as 
muons (within the accelerator acceptance) at the end of 
the cooling channel. Figure [2] shows the expected muon 
yields for the Study 2a geometry target configuration, as 
a function of the kinetic energy of the incoming proton 
beam. These yields are normalised to the total number 
of generated protons as well as the initial beam kinetic 
energy. The results show that the optimum proton beam 
energy for producing muons is approximately 8 GeV. For 
a given proton beam energy, there is significant variation 
in the calculated muon yield, depending on what version 
of the simulation codes are used, as well as on the dif- 
ferent assumptions made about the proton beam bunch 
spacing (0 or 3 ns). For the former, this can only be 
improved by updates to the simulation codes. 



B. Energy deposition 

Figure [3] shows the distribution of energy deposition in 
the Study 2a target station using the FLUKA simulation 
code, which agrees rather well with those results obtained 
using MARS [93 . We can see that the energy deposition 



in the first few superconducting coils is too high (53, 6 
and 1 kW) , meaning that they will immediately quench. 
To eliminate this problem, the shielding volume needs to 
be effectively doubled in size Previous studies have 
shown that the forces between the superconducting coils 



are very high [11 



and in the present geometry there is 
no space for support structures between individual coils, 
which needs to be addressed. Another issue is that the 
normal conducting magnets also experience very high en- 
ergy deposition rates (142 and 90 kW). This problem is 
avoided by removing them completely and replacing the 
volume with more shielding. This means that the mag- 
netic field in the beam-jet interaction region is decreased 
from 20 T down to 15 T. Also note that the energy de- 
position in the mercury pool is quite low (13 kW). This 
means that the volume of the pool reservoir needs to be 
increased in order for it to act as an effective beam dump 
for protons that do not interact with the mercury jet. In 
the next section we explore what effect these necessary 
geometry changes have on both the distribution of the 
energy deposition and on the accepted muon yield effi- 
ciencies. 



IV. INCREASED SHIELDING GEOMETRY 

Figure ID shows the new target station geometry config- 
uration that has increased shielding (approximately dou- 
bled in size) to protect the superconducting (SC) coils, 
with the normal conducting magnets removed, as well as 
a larger mercury pool reservoir (88 cm < z < 367 cm, 
y < —15 cm, r < 45 cm). Table |T] provides the geo- 
metrical and current density parameters of the SC coils. 
The SC coils are arranged in triplets, each corresponding 
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FIG. 3. Distribution of the deposited energy density (J/cc) 
per beam pulse (50 Hz repetition rate) in the Study 2a target 
system. Also shown is the estimated deposited power within 
various regions. 



The beam pipe is modelled as a 2 cm-thick stainless 
steel tapered volume, which defines the inner bore of the 
decay vacuum region. The beam pipe section just above 
the mercury pool surface is removed between z = 88 
and 200 cm to allow the mercury jet (and non-interacting 
protons) to enter the pool unhindered. The beryllium 
window, which separates the jet-beam interaction region 
from the rest of the decay channel vacuum, is moved 
further upstream to coincide with the edge of the first 
cryostat module. It consists of a 0.5 cm gap containing 
rapidly flowing He gas for cooling sandwiched between 
two 0.5 cm Be layers. The iron yoke plug has been re- 
moved to allow for space for the mercury jet nozzle in- 
jection and return flow system (which is ignored in the 
simulation) . The mercury jet is still modelled as a simple 
cylinder of radius 4 mm, tilted at approximately 100 mrad 
with respect to the magnetic z axis. 



to one cryostat module. Gaps are introduced between 
neighbouring coils to provide space for the cryostat mod- 
ules and supporting structures to hold them in place (the 
cooling components and their internal shielding are not 
included in the simulation). Each cryostat is protected 
by large shielding volumes, made from tungsten beads 
(60%) with an assumed low-grade density of 15.8 g/cc, 
and fast- flowing helium gas (40%), which must be used 
as the cooling agent, since water reacts with tungsten. 
Each shielding section, separated by 20 cm gaps, is sur- 
rounded by stainless steel container vessels (each with a 
thickness between 2 and 10 cm) that must support the 
^200 tonne weight while limiting stresses and deforma- 
tions to acceptable values p^ . 

TABLE I. Parameters defining the superconducting (SC) coils 
for the new increased shielding geometry shown in Fig. [l] zq 
is the initial z position, Az specifies the length along ri 
is the inner radius, Ar is the radial thickness, while 7 is the 
average current density. The last column specifies the coil 
materials used in the FLUKA simulation, where SCon is the 
compound Cu (54%), Nb (24%), Ti (12%) & Kapton (10%). 
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355.0 


120.0 
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FIG. 4. Schematic of the updated target geometry configura- 
tion with increased shielding. The superconducting magnets 
are labelled SCn, where n = 1 to 12. 



A. Pion and muon yields 

The removal of the normal conducting magnets im- 
plies that the peak magnetic field in the jet-beam inter- 
action region is decreased from 20 T to 15 T. To ensure 
that secondary pions from the target maintain their mag- 
netic rigidity, the inner bore radius of the beam pipe at 
this location is increased from 7.5 cm to 10 cm. Figure [5] 
shows a comparison of the axial magnetic field profile 
between the original Study 2a geometry and the new in- 
creased shielding geometry. The latter has a much wider 
magnetic field profile, which helps to capture much more 
secondary charged particles from the target. This can 
be clearly seen in Fig. [HI which shows the charged av- 
eraged yield of useful pions and muons from the target, 
normalised to the total number of generated protons and 
the initial beam kinetic energy. Using the FLUKA simu- 
lation package, these yields are calculated by finding the 
number of pions and muons (of both signs) that pass a 



4 



transverse plane 50 m downstream from the beam-jet in- 
teraction region, within the decay channel aperture that 
has a bore radius of 30 cm. In order to obtain the approx- 
imate number of useful muons for the Neutrino Factory 
that pass through the cooling channel, we require that 
these particles have kinetic energies in the range of 40 to 
180 MeV. 

In addition to the baseline mercury jet case, the 
yields from alternative target materials are also shown 
in Fig. ini using the simulation parameters outlined in 
Sec. |TT1 At intermediate and high beam kinetic energies 
(above 5 GeV), the mercury jet provides the best yields, 
although the performance of the solid tungsten target 
matches this very closely. The results show that a rel- 
ative gain of 10% in the particle yield is possible if the 
proton beam energy is increased to 15 GeV. Note that 
we do not see a dramatic reduction in the yield for the 
tungsten powder jet, which is assumed to have an effec- 
tive density of 50% of solid tungsten. The reason for this 
is that, even though less protons will interact with the 
powder target, the amount of re- absorption of secondary 
particles inside the target will also be lower, giving an 
overall figure of merit comparable to either the mercury 
jet or solid tungsten case. At low beam kinetic energy 
(below 5 GeV) , there is an indication that the liquid gal- 
lium target gives the best yield compared to the other 
materials. 
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FIG. 5. Axial magnetic field profile distributions for the Study 
2a and new increased shielding geometries. 
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FIG. 6. The charged-averaged accepted pion and muon yield 
per proton per GeV for various targets in the new increased 
shielding geometry. Also shown are the equivalent yields for 
the mercury jet target in the Study 2a geometry. 



duced the power deposition in the SC coils. Table [TT] 
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FIG. 7. Distribution of the deposited energy density in the 
increased shielding geometry configuration for the mercury jet 
target. Also shown is the estimated power deposition within 
various regions (blue italics denote the power deposited in the 
beam pipe). 



B. Energy deposition 

We have seen that the change to the target station ge- 
ometry has had a positive effect to the useful muon yield. 
The situation regarding the energy deposition for the su- 
perconducting coils is also improved. Figure [7] shows the 
radiation dose rate in the new target station geometry, 
obtained using the FLUKA simulation code, which shows 
that the increase to the shielding has dramatically re- 



provides a detailed breakdown of the deposited power 
in various sub-regions of the new geometry, where the 
uncertainties are estimated by using different initial ran- 
dom number seeds for the simulation. The first SC coil 
experiences a total power of 0.4 kW, with the other coils 
experiencing much lower values, giving a total approxi- 
mately equal to 0.6 kW. These correspond to energy den- 
sities below O.lmW/g, which permits at least a 10 year 
(of 2 X 10^ s each) operational lifetime against radiation 
damage to their organic (MgO) insulators, assuming the 
maximum permitted radiation dose is lOMGy [l3| . 
About 2.4 MW is deposited in the combined shielding 
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TABLE II. Power deposition in various regions of the increased shielding geometry configuration for the mercury jet, as well 
as for alternative target materials. For the solid tungsten bars, the pool reservoir is replaced by more shielding. 

Region Deposited power (kW) 







Hgjet 


Ga jet 


W powder jet 


W solid 


SC coils 1-12 




0.57 ±0.05 


0.67 ±0.06 


0.62 ± 0.06 


0.55 ±0.06 


Lower shielding SC 1-3 (r < 50 cm, 


z < 83 cm) 


1284.4 ± 8.3 


1034.9 ±8.1 


1154.3 ±9.0 


1282.6 ± 7.1 


Lower shielding SC 1-3 (r < 50 cm, 


z > S3 cm) 


234.2 ± 3.5 


318.5 ±3.6 


284.8 ±4.0 


348.2 ±6.2 


Upper shielding SC 1-3 (r > 50 cm) 




58.3 ±0.8 


82.6 ± 1.1 


75.8 ± 1.0 


41.9 ±0.5 


Shielding for SC 4-6 




38.0 ±1.7 


45.2 ± 2.0 


38.8 ± 1.9 


26.1 ± 1.2 


Shielding for SC 7-9 




11.0 ±0.7 


11.9 ±0.8 


10.8 ±0.8 


8.1 ±0.6 


Shielding for SC 10-12 




7.4 ±0.7 


8.1 ±0.7 


7.1 ±0.8 


5.0 ±0.5 


Beam pipe up to 2 = cm 




352.3 ± 2.9 


230.8 ± 2.3 


303.2 ± 3.3 


303.0 ± 1.9 


Beam pipe from z = cm to end of taper 


397.6 ± 3.6 


499.1 ±4.3 


428.7 ±3.8 


338.8 ±4.1 


Beam pipe from end of taper 




21.7 ±0.9 


24.0 ± 1.0 


21.1 ± 1.0 


14.6 ±0.8 


Lower shielding vessel for SC 1-3 (r 


< 50 cm) 


7.6 ±0.2 


12.5 ± 0.3 


9.7 ±0.2 


5.7 ±0.2 


Upper shielding vessel for SC 1-3 (r 


> 50 cm) 


6.0 ±0.1 


9.3 ±0.2 


7.8 ±0.2 


4.3 ±0.1 


Shielding vessel for SC 4-6 




3.5 ±0.3 


4.5 ±0.3 


3.6 ±0.3 


2.4 ±0.2 


Shielding vessel for SC 7-9 




0.8 ±0.1 


0.8 ±0.1 


0.7 ±0.1 


0.6 ±0.1 


Shielding vessel for SC 10-12 




0.5 ±0.1 


0.6 ±0.1 


0.5 ±0.1 


0.3 ±0.1 


Pool reservoir container 




10.5 ±0.3 


17.1 ± 0.4 


14.0 ± 0.4 




Pool reservoir 




460.8 ±9.7 


814.1 ± 10.6 


655.3 ± 11.5 




Jet/Target 




416.8 ±2.4 


167.3 ± 1.0 


298.7 ±2.3 


1018.5 ± 5.2 


Be window 




8.9 ±0.1 


6.3 ±0.1 


8.4 ±0.1 


5.1 ±0.1 


Total 




3320.7 ± 14.4 


3288.1 ± 15.0 


3323.8 ± 16.3 


3405.6 ± 11.8 



and beam pipe sections. Removing the total heat load 
from these materials will prove quite challenging. The 
target and mercury pool receives a total power deposi- 
tion approximately equal to 0.9 MW, which must be dis- 
sipated in a heat exchanger in the mercury flow return 
loop. The downstream Be window receives a large energy 
deposition of 9 kW. About 250 kW will continue into the 
downstream pion-muon transport system. This will con- 
tain scattered, high energy protons that must be removed 
by a chicane and absorber system in order to s top radi- 
ation damage to other accelerator components jl4| . 

Table [TT] also shows the deposited power in the vari- 
ous regions when alternative materials are used for the 
target. In general, most of the differences are seen in 
the energy deposition in the target itself, the shielding 
and the liquid/powder collection pool. Less protons in- 
teract with the liquid gallium jet target, going directly 
into the pool reservoir. This very slightly reduces the 
overall power deposited in the shielding and beam pipe 
sections. The reduced density in the tungsten powder jet 
also means that more energy is deposited in the powder 
collection reservoir (assumed to have the same geometry 
as the mercury pool), although not as much as was the 
case for the gallium target. In contrast, the solid tungsten 
target itself experiences a much larger energy deposition, 
corresponding to a quarter of the total beam power. Ex- 
perimental work using fast, high current pulses passing 
through tungsten wires has demonstrated that a 20 cm- 



long, 2 cm diameter tungsten target will be able to with- 
stand the peak stresses from such an energy deposition, 
allowing a target lifetime of at least three years hj. The 
removal of the reservoir pool for the solid target means 
that about 5% more energy is deposited in the shielding, 
with a slight reduction in the radiation dose absorbed by 
the downstream beam pipe. 



V. FAILURE MODES 

It is important to know how the 4MW total beam 
power will be distributed within the new target station 
geometry for specific (baseline) failure modes: when the 
magnetic field fails, when the mercury jet stops flowing, 
and when both happen together. Table IIIII shows com- 
parisons of the average power deposition in various re- 
gions of the target station between normal and failure- 
mode operating scenarios. 

Under normal conditions, non-interacting protons from 
the beam have a trajectory that enters the mercury pool 
reservoir. This is illustrated by the case when there is 
no mercury jet target, in which more than half of the 
total beam power (2.6 MW) is deposited in the pool. 
This will produce significant agitation of the pool sur- 
face with splashes of radial velocities expected to ap- 
proach 50 ms""'^ [l6j . In contrast, when there is no mag- 
netic field present to steer the proton beam (and any 
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TABLE III. Power deposition in various regions of the increased shielding geometry configuration for the mercury jet for a 
range of operational failure modes. 

Region Deposited power (kW) 







Nominal 


No Hg jet 


No B 


No B & No Hg jet 


SC coils 1-12 




0.57 ±0.05 


0.43 ± 0.05 


0.50 ± 0.05 


0.51 ±0.05 


Lower shielding SC 1-3 (r < 50 cm, 


z < 83 cm) 


1284.4 ± 8.3 


42.0 ±0.3 


842.8 ± 5.8 


842.9 ±5.5 


Lower shielding SC 1-3 (r < 50 cm, 


2 > 83 cm) 


234.2 ±3.5 


668.3 ±2.9 


1370.3 ± 9.4 


1371.2 ±9.4 


Upper shielding SC 1-3 (r > 50 cm) 




58.3 ± 0.8 


120.7 ±0.8 


61.6 ±0.8 


61.7 ±0.7 


Shielding for SC 4-6 




38.0 ± 1.7 


15.0 ±1.0 


23.8 ± 1.4 


23.9 ± 1.3 


Shielding for SC 7-9 




11.0 ±0.7 


0.8 ±0.3 


3.1 ±0.5 


3.0 ±0.4 


Shielding for SC 10-12 




7.4 ±0.7 


0.2 ±0.1 


1.6 ±0.3 


1.6 ± 0.4 


Beam pipe up to 2 = cm 




352.3 ± 2.9 


0.2 ±0.1 


0.7 ±0.1 


0.7 ±0.1 


Beam pipe from z = cm to end of taper 


397.6 ± 3.6 


218.8 ± 1.7 


1013.5 ±5.0 


1012.8 ± 5.5 


Beam pipe from end of taper 




21.7 ±0.9 


0.5 ±0.2 


4.6 ±0.4 


4.6 ±0.5 


Lower shielding vessel for SC 1-3 (r 


< 50 cm) 


7.6 ±0.2 


12.0 ±0.2 


8.0 ±0.2 


8.0 ±0.2 


Upper shielding vessel for SC 1-3 (r 


> 50 cm) 


6.0 ±0.1 


11.8 ±0.1 


6.3 ±0.1 


6.2 ± 0.1 


Shielding vessel for SC 4-6 




3.5 ±0.3 


1.9 ±0.2 


2.6 ±0.3 


2.6 ± 0.2 


Shielding vessel for SC 7-9 




0.8 ±0.1 


0.1 ±0.1 


0.2 ±0.1 


0.2 ± 0.1 


Shielding vessel for SC 10-12 




0.5 ±0.1 


< 0.1 


0.1 ±0.1 


0.1 ±0.1 


Pool reservoir container 




10.5 ±0.3 


10.8 ±0.1 


14.1 ±0.2 


14.2 ± 0.2 


Pool reservoir 




460.8 ±9.7 


2603.2 ±4.8 


235.9 ± 2.5 


237.5 ±2.4 


Jet/Target 




416.8 ± 2.4 




1.9 ±0.1 




Be window 




8.9 ±0.1 


0.1 ±0.1 


0.3 ±0.1 


0.3 ±0.1 


Total 




3320.7 ± 14.4 


3707.0 ± 6.0 


3591.9 ± 12.5 


3591.9 ± 12.6 



secondary charged particles), the energy deposition for 
the combined mercury jet and pool system decreases by 
roughly a factor of three. The proton beam almost com- 
pletely misses the mercury jet target and hits only part of 
the mercury pool, instead dramatically increasing the en- 
ergy deposition in the surrounding shielding (1.6 MW to 
2.3 MW) and the nearby beam pipe section (from 0.4 MW 
to 1 MW) . For the scenario when there is no mercury jet 
and no magnetic field, the results are essentially identi- 
cal to the case when only the magnetic field is turned 
off, owing to the fact that the proton beam will miss the 
target in both cases. 



VI. RADIATION SAFETY 

Any construction of the target station must take into 
account the safety requirements of the surrounding en- 
vironment to prevent radiation contamination of the soil 
(and ground water). The entire target station must be 
enclosed within a concrete shielding structure of adequate 
thickness to stop radiation from escaping; the effective 
total radiation dose must not exceed 1 mSv per year [l7j , 
which is equivalent to a continual dose of 0.1 /zSv/hr. 

Figure |8] shows the estimated total radiation ambient 
dose equivalent from the 4 MW target station after a total 
irradiation time of 2 x 10^ s (1 year), using the FLUKA 
simulation code and AMB74 conversion factors [l^. The 



target station is surrounded by a concrete shielding struc- 
ture (tunnel) that is modelled as three connecting cylin- 
drical sections, extending from z = —8 m up to 2: = 30 m, 
with an outer radius of 6 m. This tunnel creates a barrier 
between the surrounding rock, assumed to be molasse 
soil, and the radiation generated by the target system. 
The chemical composition (with mass fractions in paren- 
thesis) assumed for the concrete is O (51.1%), Si (35.8%), 
Ca (8.6%), Al (2.0%), Fe (1.2%), H (0.6%), C (0.4%) and 
Na (0.3%), while the composition of the molasse soil is 
taken to be O (49.2%), Si (19.8%), Ca (9.3%), Al (6.4%), 
C (4.9%), Fe (4.1%), Mg (3.5%), K (1.9%), Na (0.6%), 
Mn (0.2%) and Ti (0.1%) The densities assumed 

for the concrete and molasse are 2.35 g/cc and 2.4g/cc, 
respectively. For a range of cooling decay times, it can 
be seen that no radiation escapes the concrete shielding, 
which means that there will be minimal activation of the 
surrounding soil and ground water. We can also infer 
that remote handling will be mandatory for maintaining 
the target system; even at a radial distance of 2 m from 
the interaction region, the activity after 1 year of cooling 
is of the order of lOmSv/hr, which greatly exceeds the 
safety limit for radiation workers (2GmSv/yr) [l7| . 
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d) 1 year cooling 
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FIG. 8. Ambient dose equivalent rates for the target station, concrete shielding tunnel and surrounding underground rock. 



VII. SUMMARY 



We have presented FLUKA and MARS simulation 
studies of the pion production and energy deposition in 
the Neutrino Factory target station. Compared to the 
original Study 2a geometry, a doubling of the shielding 
volume is needed to protect the superconducting coils 
from too much radiation damage. In addition, a reduc- 
tion in the focusing magnetic field from 20 T down to 15 T 
does not affect the pion production efficiency, provided 
the inner bore radius is increased. Alternative target ma- 
terials were also investigated, such as liquid gallium, and 
powdered and solid tungsten. Each of these offer com- 
parable muon yields, with some differences observed in 
the power deposited in the target, collection pool and 
surrounding shielding. We have investigated what effect 



various operational failure modes have on the deposited 
power in the target station. Finally, we have provided 
estimates of the amount of concrete shielding that will 
be needed to protect the environment from the high ra- 
diation generated by the target station, with remote han- 
dling mandatory for any maintenance work. 
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